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SECTION  I 


INTRODUCTION 

The  objective  of  this  contract  research  was  to  provide  more  insight  into  die  sensitivity 
of  nitroaromatic  compounds  toward  shock  and  impact  by  means  of  ab  initio  molecular  orbital 
calculations  on  nitro  and  oci-nitro  compounds.  A  specific  task  was  to  explore  the  relevance 
of  oci-nitro  derivatives  in  the  decomposition  pathway  of  nitroaromatics.  The  target  systems 
selected  by  the  Wright  Laboratory  /  Armament  Directorate  of  were  trinitrotoluene  (TNT), 
trinitrophenol  (picric  acid),  and  trinitroaniline  (picramide)  with  the  understanding  that 
computational  rigor  might  be  limited  for  these  large  systems.  The  basis  for  the  research  was 
a  tantomerism  study  on  o-nitrophenol  reported  by  Politzer  and  coworkers  (Reference  1).  Our 
approach  to  the  relevance  of  nitro  aci-  nitro  tautomerism  has  been  to  delineate  the 

computational  aspect  for  this  process  "from  the  ground  up."  Therefore,  we  have  studied  the 
parent  tautomeric  processes  (nitro  **  oci-nitro,  imine  «*  enamine,  and  keto  <•*  enol)  in  systems 
as  simple  as  nitromethane,  in  composite  nitroethylene  systems,  and  in  substituted 
nitroaromatics.  This  approach  enabled  us  to  set  computational  error  bars  for  structural  and 
energetic  data  and  to  make  comparisons  with  experimental  data.  In  addition  to  investigating 
the  relevance  of  oci-nitro  compounds,  the  primary  focus,  we  have  also  touched  upon  their 
involvement  in  reaction  pathways.  Calculations  have  been  conducted  on  a  large  number  of 
systems  that  collectively  lead  to  important  observations.  For  clarity,  we  discuss  in  separate 
sections  the  parent  tautomerisms,  composite  systems,  and  nitroaromatics.  In  this  section  we 
briefly  describe  relevant  background  data. 

1  BACKGROUND 

The  question  this  project  addresses  is:  why  is  picric  acid  mote  sensitive  and  picramide 
less  sensitive  than  TNT,  as  measured  by  impact  drop  height,  (Le.  87  cm  for  picric  acid,  177 
cm  for  picramide,  and  160  cm  for  TNT  (Reference  2))?  Thermochemical  decompositions  of 
TNT  by  isothermal  differential  scanning  calorimetric  (DSC)  analyses  gave  primary  kinetic 
deuterium  isotope  effects  (KDIE)  of  1.66  (24S-269  °C)  at  induction  and  of  1.35  during  the 
decay  (References  3  and  4),  which  suggest  that  C — H  bond  rupture  occurs  in  the  rate¬ 
determining  step.  The  induction  activation  energy  for  the  decomposition  of  the  liquid  was 
estimated  (from  Arrhenius  parameters)  at  46.5  kcal/mol  with  a  decay  activation  energy  of  29 
kcal/mol  (Reference  4).  C — NOj  Bond  fission  is  the  preferred  gas-phase  decomposition 
process  of  mono-  and  dinitro toluenes,  using  laser-powered  homogeneous  pyrolysis,  with 
activation  energies  (from  Arrhenius  parameters)  of  ca.  70  kcal/mol  (Refererence  5).  Incipient 
formation  of  anthranil  occurs  in  the  pyrolysis  of  o-nitro  toluene,  using  a  single  shock  pulse 
tube  (Reference  6);  2,4-dinitroanthranil  is  a  known  TNT  decomposition  product  (Reference 
7).  Picosecond  transient  absorption  studies  of  2,4-  and  2,6-dinitrotoluene  solutions  have 
suggested  the  formation  of  the  ad-forms  (Reference  8).  Electron  paramagnetic  resonance 
(EPP '  studies  on  the  thermochemical  decomposition  of  liquid  TNT  showed  the  formation  of 
aryi  'enzyloxynitroxide  radicals  by  intermolecular  coupling  (Reference  9).  Using  EPR, 
activation  energies  of  40.9  (30.2)  kcal/mol  for  the  induction  (decay)  were  determined  in  the 
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thermochemical  decomposition  of  TNT  (Reference  10),  which  are  similar  to  those  determined 
by  DSC 

Nitro  **  ad-nitro  tautomerism  was  used  by  Politzer  and  coworkers  (Reference  1)  in  a 
theoretical  (HF/3-21G)  study  to  interpret  differences  in  sensitivities  between  hydroxy  and 
amino  substituted  o-nitrobenzene.  Cox  and  Hillier  (Reference  11)  reported  on  the  aci-form  of 
TNT  at  a  rudimentary  ab  initio  level  (HF/STO-3G).  Calculated  geometries  and  energies  of 
nitro  compounds  are,  however,  very  sensitive  to  the  theoretical  levels  employed.  For 
example,  Head-Gordon  and  Pople  (Reference  12)  calculated  a  4.55  (7.2)  kcal/mol  N02- 
rotation  barrier  at  MP2/6-311G**  (HF)  (Reference  13)  for  nitrobenzene  while  the 
experimental  value  is  2.9  kcal/mol;  the  HF-MP2/6-31G*  geometry  difference  is  small 
(Reference  14).  This  barrier  amounts  to  4.83  kcal/mol  at  QCISD0D/6-311G**  for 
nL  oethylene  at  which  level  it  is  essentially  the  same  as  the  experimental  value  of  5.05 
kcal/mol;  the  HF-MP2/6-31G*  geometry  difference  amounts  to  5%  (Reference  15).  Bond 
dissociation  energies  are  even  more  difficult  to  calculate.  For  example,  the  MRCI/6- 
31G*//CAS/6-31G*  calculated  CHj-NOj  -4  CHs-ONO  rearrangement  barrier  is  10  kcal/mol 
too  high  (Reference  16).  These  examples  illustrate  the  need  for  careful  calibration  of  the 
theoretical  methods. 

2  COMPUTATIONAL  METHODS 

Ab  initio  molecular  orbital  calculations  were  carried  out  using  the  GAUSSIAN  90/92 
suite  of  programs  (Reference  17).  Structures  were  fully  optimized  within  given  symmetries 
with  the  split  valence  HF/3-21G  basis,  the  heavy  atom  d-polarized  HF/6-31G*  basis  (nitro- 
aromatics)  and  MP2(full)/6-31G*  (all  other  species),  which  includes  the  effects  of  all  electron 
correlation  using  second-order  M0ller-Plesset  perturbation  theory.  The  identity  of  each 
structure  was  established  at  its  highest  level  by  determining  the  number  of  negative  eigen¬ 
values  of  the  force  constant  matrix.  These  structures  are  also  the  basis  for  the  frozen  core 
MP4(fc)  single  point  higher  level  calculations  (for  the  smallest  systems)  that  employ  up  to 
essentially  a  triply  split  basis  set  (6-31 1G)  with  d-  and  hydrogen  p-polarization  functions  (**) 
and  diffuse  p-  and  s-functions  (++);  MP4(fc)/6-31G*  was  used  for  the  nitroethylenes  and 
MP2/6-31G*  for  the  nitroaroraatics.  For  the  smallest  species  Gaussian  1  (Gl)  energies 
(equation  1)  are  reported  except  that  0.95  scaled  MP2/6-31G*  ZPEs  are  used. 

Gl  =  MP4/6-311G**  +  AE(+)  +  AE(2df)  +  AE(QCI)  +  AE(HCL)  +  scaled  ZPE  [1] 

Bonding  properties  were  investigated  with  Bader's  topological  one-electron  density 
analysis  (Reference  18)  using  optimized  HF/6-31G*  or  MP2(full)/6-31G*  wave  functions. 

The  one-electron  density  distributions  p(r)  were  analyzed  with  the  aid  of  the  gradient  vector 
field  Vp(r)  and  the  Laplacian  V2p(r),  which  also  determines  the  regions  in  space  wherein 
electronic  charge  is  concentrated  or  depleted.  Bond  critical  points  have  a  minimum  value  in 
p(r)  along  the  maximum  electron  density  path  connecting  two  nuclei  and  are  maxima  in  all 
other  directions.  The  ellipticity  e  at  such  a  critical  point  describes  the  spatial  symmetry  of  a 
bond,  whereas  the  total  energy  density  indicates  a  bond  to  be  covalent  (H(r)  <  0)  or  ionic 
(H( r)  >  0).  The  electron  population  of  elements  in  structures  (charges)  were  obtained  by 
integration  of  the  charge  density  within  their  respective  basins  as  defined  by  the  zero-flux 
surfaces. 
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SECTION  n 


NITRO  ^  aci-NTTRO,  KETO  ENOL,  AND  IMINE  EN AMINE  TAUTOMERISM 

1.  INTRODUCTION 

In  order  to  establish  the  validity  of  ori-nitro  compounds  in  the  context  of  this  research, 
it  is  of  importance  to  determine  at  which  levels  of  theory  their  structures  and  energies  are 
calculated  accurately.  Nitronic  acids  are  well  known  in  organic  chemistry,  but  a  large  energy 
difference  of  21.8  kcal/mol  for  the  nitromethane  and  oct-nitromethane  tautomers  was  reported 
at  MP2/6-31G*  (Reference  19),  which  is  for  most  organic  systems  a  respectable  level  of 
theory.  Similar  observations  are  made  for  keto  -  enol  tautomers,  but  here  it  has  been  shown 
that  higher  levels  of  theory  reduce  the  energy  difference  between  acetaldehyde  and  vinylalco* 
hoi  from  16.4  kcal/mol  at  MP2/6-31G*  to  10.4  kcal/mol  at  MP3/6-31++G**  (Reference  20). 
We  chose  to  use  G1  theory  to  set  a  uniform  standard  and  to  determine  the  dependence  of 
structures  and  energies  on  the  levels  of  theory  for  the  tautomerisms  shown  below. 


nrtro/acHiitro 


keto/enol 


imine/enamine 
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2. 


NITRO  ori-NITRO  TAUTOMERISM 


We  consider  the  nitromethane  **  aci-nitromethane  tautoraerism.  Their  structures  and 
that  of  their  common  anion  were  optimized  at  all  levels  of  theory  of  which  the  HF  and 
MP2/6-31G*  parameters  are  displayed  in  Figure  1.  Total,  relative,  and  ionization  energies  are 
listed  in  Tables  1-3,  respectively.  The  properties  of  the  bond  critical  points  of  some  structures 
are  provided.  Atomic  properties  are  shown  in  Figure  2. 

a.  Discussion 

The  MP2/6-31G*  geometry  of  nitromethane  is  in  good  agreement  with  those  obtained  from 
microwave  spectroscopy  (assigned  to  2,  Reference  21)  and  single  neutron  diffraction  data  at 
IS  K  (assigned  to  1,  Reference  22),  but  deviates  slightly  from  a  single  crystal  X-ray  structure 
at  228  K  (Reference  23),  which  is  just  below  the  melting  point  These  data  do  suggest  that 
the  N — O  bond  lengths  calculated  at  MP2/6-31G*  are  ca.  0.015  A  too  long.  MP2  has  a 
tendency  to  exaggerate  the  loss  of  charge  density  in  the  bonding  region,  and  this  effect  is 
strongest  for  the  NO  bonds,  i.e.,  the  MP2-HF/6-31G*  differences  in  bond  lengths  are  AdNO  = 
0.051  A.  Similar  observations  are  made  for  the  harmonic  frequencies,  which  underscores  that 
the  N02-group  is  the  least  well  described  at  MP2/6-31G*.  Electron  density  analysis  of  the 
bonds  suggests  that  there  is  little  if  any  hyperconjugation  between  the  CH3  and  N02  groups, 
which  is  in  line  with  the  absence  of  a  rotation  barrier.  From  the  G1  energies  an  atomization 
energy  (ZDJ  of  -  570.69  kcal/mol  is  calculated  for  nitromethane,  which  compares  very  well 
with  -  573.13  kcal/mol  that  is  deduced  from  experimental  heats  of  formation  (Reference  24). 

The  cis  and  trans  forms  of  aci-nitromethane  differ  in  the  orientation  of  the  OH  relative 
to  the  NO  unit  The  trans  conformer  has  the  higher  dipole  moment  and  is  the  transition 
structure  for  N — OH  bond  rotation  (at  MP2/6-31G*)  with  a  barrier  of  6.8  kcal/mol  at  Gl. 

The  MP2-HF  geometrical  differences  are  largely  in  the  C=N  and  N — OH  bonds.  The  electron 
density  analysis  shows  aci-nitromethane  to  be  strongly  polarized  with  an  essentially  neutral 
nitrogen  atom.  The  energy  difference  between  this  tautomer  and  nitromethane  amounts  to 
14.1  kcal/mol  at  Gl.  This  energy  difference  is  significantly  larger  at  lower  levels  of  theory. 
For  example,  AE  is  19.1  kcal/mol  at  MP4/6-311+4G**  +  scaled  ZPE,  22.9  kcal/mol  at 
MP2/6-31G*  (ZPE  uncorrected),  20.3  kcal/mol  at  HF/6-31G*,  and  0.3  kcal/mol  at  HF/3-21G. 

Another  measure  for  the  accuracy  of  the  calculations  is  the  determination  of  proton 
affinities.  The  Gl  calculated  proton  affinity  of  the  nitromethide  anion  is  355.2  kcal/mol 
which  compares  very  well  with  the  experimentally  determined  ionization  energy  of  356.4 
kcal/mol  (gas  phase,  pulsed  ICR,  Reference  25)  for  nitromethane.  The  electron  density 
analysis  supports  the  contention  that  the  strongly  negatively  charged  oxygens  are  the  most 
receptive  toward  protonation,  which  would  render  aci-nitromethane.  The  anion  is  a  potential 
Y-aromatic  system.  However,  the  methylene  group  pyramidalizes  (15.7°)  using  MP2/6-ll-tG* 
at  which  level  the  inversion  barrier  is  a  small  32  cal/mol. 
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Optimized  Nitromethane  and  acj-Nitromethane  Tautomers  and  their  Common  Anion.  The  MP2/6-31G*  Geometrical 
rs  are  underlined.  The  Lower  Values  are  Those  at  HF76-31G*.  The  Top  Values  in  Italics  are  Those  at  MP2/5-3 1 1+G*. 
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Figure  2.  Atomic  and  Group  Charges  for  die  Nitromethane  and  act-Nftromethane  Tautomers  and  their 


3. 


KETO  **  ENOL  AND  IMINE  ENAMINE  TAUTOMERISM 


The  parent  acetaldehyde  vinyl  alcohol  and  acetaldimine  vinyl  amine  tautomerisms 
are  considered.  Their  structures  and  those  of  their  conformers  and  anions  were  optimized  at 
all  levels  of  theory  of  which  the  HF  and  MP2/6-31G*  parameters  are  displayed  in  Figure  3. 
Total*  relative,  and  ionization  energies  are  listed  in  Tables  4-6,  respectively.  Atomic 
properties  are  shown  in  Figure  4.  Harmonic  frequencies  are  listed  in  Table  7. 

a.  Discussion 

The  MP2/6-31G*  structures  compare  well  with  those  available  experimentally.  This  is 
the  case  for  acetaldehyde  (Reference  26),  vinyl  alcohol  (microwave  spectrum.  Reference  27), 
acetaldimine  (microwave  spectrum  (Reference  28)  and  gas  phase  (Reference  29)  and  matrix 
1R  spectra  (Reference  30))  and  vinyl  amine  (microwave  spectrum.  Reference  31).  The  HF- 
MP2/6-31G*  difference  in  geometrical  parameters  is  largest  for  C — O,  0=0,  C — N,  C=N,  and 
C=C  bonds.  The  rms  differences  between  experimental  and  0.95  scaled  MP2/6-31G* 
harmonic  frequencies  are  28  cm'1  for  acetaldehyde  (Reference  32),  27  cm'1  for  vinyl  alcohol 
(Reference  33),  30  cm'1  for  acetaldimine  (Reference  29),  52  cm*1  for  vinyl  amine  (Reference 
34),  which  may  suggest  that  acetaldimine  and  vinyl  amine  are  the  least  well  described  experi¬ 
mentally. 

Vinyl  alcohol  has  two  conformers  of  which  the  syn  form  is  0.9  kcal/mol  (Gl)  more 
stable  than  the  anti  conformer.  In  contrast,  ann-acetaldimine  is  0.9  kcal/mol  more  stable  than 
the  syn  form  (MP4/6-311+4G**),  which  agrees  with  the  experimental  value  of  ca.  2  kcal/mol 
(Reference  35),  and  corrects  earlier  theoretical  data  (Reference  36).  Vinyl  amine  is  calculated 
to  be  non-planar  with  an  N-inversion  barrier  of  1.4  kcal/mol  (MP4/6-311++G**),  which 
agrees  well  with  the  best  experimental  estimates  of  1.1  kcal/mol  (Reference  37).  The  anion 
of  acetaldimine  prefers  the  syn  conformation  over  the  trans  form  by  0.8  kcal/mol  (Gl). 

The  Gl  energy  difference  between  the  acetaldehyde  and  vinyl  alcohol  tautomers  is 
10.9  kcal/mol,  which  compares  well  with  the  estimated  gas  phase  value  of  9.9  kcal/mol.  A 
much  smaller  energy  difference  of  3.8  kcal/mol  (Gl)  is  calculated  for  the  acetaldimine  and 
vinyl  amine  tautomers,  for  which  no  experimental  value  is  available. 

The  theoretically  estimated  atomization  energy  of  -  641.4  kcal/mol  (Gl)  for  acetalde¬ 
hyde  is  in  good  agreement  with  the  experimental  value  of  -  642.5  kcal/mol  (Reference  38). 

The  Gl  theoretical  estimate  for  the  atomization  energy  of  acetaldimine  is  -  694.3  kcal/mol, 
which  differs  significantly  from  the  experimental  value  of  -  708.7  kcal/mol.  The  theoretical 
estimate  for  the  ionization  of  acetaldehyde  is  365.0  kcal/mol,  which  is  in  good  agreement  with 
the  gas  phase  value  of  366.4±3  kcal/mol  (Reference  39).  Likewise,  the  calculated  ionization 
energy  of  354.1  kcal/mol  for  syn- vinyl  alcohol  compares  well  with  the  experimental  value  of 
356.6  kcal/mol.  The  Gl  estimated  ionization  energies  are  377.3  kcal/mol  for  anti - 
acetaldimine  and  373.5  kcal/mol  for  vinyl  amine,  but  no  experimental  data  are  available.  For 
comparison,  the  gas  phase  ionization  energy  of  propene  is  390±3  kcal/mol. 

Possible  internal  Coulombic  (-  +- )  stabilization  of  the  acetaldehyde  and  acetaldimine 
anions  was  evaluated  by  means  of  MP2/6-31G*  atomic  charges  which  are  displayed  in  Figure 
4.  The  negative  charge  is  similarly  distributed  over  the  CH2  and  CHNH  or  CHO  groups.  The 
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accumulation  of  charge  for  die  CHj  group,  as  based  on  the  same  group  in  vinyl 
alcohol/amine,  is  comparable  to  the  accumulation  of  charge  in  the  CHNH/CHO  group,  as 
baaed  on  the  same  group  in  acetaldehyde/acetaldimine.  From  this,  the  similar  heavy  atom 
polarization  in  die  neutrals  and  anions,  and  die  fact  that  nearly  die  full  unit  of  negative  charge 
can  be  accounted  for  in  the  hydrogens,  it  seems  that  the  anions  do  not  experience  a  strong 
charge  (-  +- )  effect  (Reference  20  and  40).  The  calculations  are  more  suggestive  of  extended 
conjugation. 
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be  MP2/6-31G*  Geometrical  Paramelan  are  Undwtinod  The  Lower  Values  are  Hwe  at  HFAS-310*. 


TABLE  4.  ABSOLUTE  ENERGIES  (in  -  hartrccs)  OP  THE  KEltVEWOL  AND  IMINE/EN AMINE  TAUTOMERS  AND  THEIR  ANIONS.1 
Compound  HF/3-21G  HF/6-310*  MP2/6-310*  MP4/6-31G*  MP4/6-311+4G*  OIk 
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TABLE  6.  IONIZATION  ENERGIES  (in  kcal/mol)  OF  THE  KETCMBNOL  Al 
Compd  HF/3-210  HF/6-31G*  MP2/6-310*  MP4/6-31G* 
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TABLE  7.  CALCULATED  MP2/6-31G*  AND  EXPERIMENTAL  FREQUENCIES  (in  cm'1)  FOR  THE 


SECTION  m 


TAUTOMERISM  IN  SUBSTITUTED  NITROETHYLENES 

1.  INTRODUCTION 

We  consider  the  ds  isomers  of  2-nitrovinyl  alcohol  2-nitro vinyl  amine  and  1-nitro- 
propene  and  their  tautomers.  Two  sets  of  tautomers  are  possible  for  each.  Important  for  this 
research  are  the  nitro  aci-nitro  tautomers  where  the  hydrogen  is  transfered  from  the  CHI, 
NHj,  and  CH,  substituents.  These  represent  symmetry  allowed  [1,5]-H  shifts  and  are 
examples  of  conjugate  tautomerism.  The  otter  tautomeric  processes  display  [131-H  transfers 
(as  in  Section  II),  and  these  are  better  conducted  in  an  acid/base  exchange  environment. 

The  structures  of  all  tautomers  were  optimized  at  all  levels  of  theory  of  which  the  HF 
and  MP2/6-31G*  parameters  are  shown  in  Figure  3  to  display  the  effects  of  basis  sets  and 
electron  correlation.  The  MP2/6-31G*  structures  for  the  tautomeric  pairs  and  their  transition 
structures  are  shown  in  Figures  6-8.  Total  and  relative  energies  are  listed  in  Tables  8  and  9, 
respectively,  and  visualized  in  Figure  9.  A  comparison  of  MP4/6-31G*  energies  between  the 
conjugated  tautomers  and  their  parents  is  given  in  Figure  10.  Atomic  properties  are  shown  in 
Figures  11  and  12,  which  display  bond  critical  point  data  and  total  electron  densities, 
respectively. 

2.  DISCUSSION 

The  structural  parameters  highlight  significant  differences  between  HF/3-21G  and 
HF/6-31G  structures,  underscore  the  importance  of  the  effects  of  electron  correlation  in  the 
geometry  optimizations,  and  compare  well  with  the  limited  literature  data  (Reference  41). 

The  cis  OH  and  NH,  substituted  nitroethylenes  are  well  set  up  for  tautomerism.  In 
fact,  the  separation  between  the  OH  hydrogen  and  the  neighboring  NO,  oxygen  in  16  is  only 
1.773  A,  and  that  between  one  of  the  NH,  hydrogens  and  the  NO,  oxygen  in  21  is  1.992  A. 
Both  the  (0)N0..JH0  (Figure  12)  and  (0)N0..HNH  interactions  represent  hydrogen  bonding 
as  confirmed  by  the  Bader  electron  density  analysis  which  identifies  their  bond  critical  points 
(Figure  11).  This  is  also  the  case  for  their  nitronic  acids  which  have  very  short  hydrogen 
bonds  of  1.616  A  for  (0)N0H...0  and  1.606  A  for  (0)N0H...NH.  1-Nitropropene  has  a 
bisected  conformation  of  the  methyl  group  (Figure  8).  The  eclipsed  form  26n  represents  the 
transition  for  CH, -rotation  (except  at  3-2 1G),  the  barrier  of  which  amounts  0.6  kcal/mol 
(MP4/6-31G*),  and  consequently  the  two  (0)N0..H,CH  distances  are  longer.  Its  nitronic 
add  tautomer  27  has  a  conformation  in  which  the  od-nitro  hydrogen  is  anti  to  the  methylene 
group  and  thereby  contrasts  the  keto  (17)  and  inline  (22)  analogues,  which  have  corresponding 
syn  conformations  for  their  aci-nitro  hydrogens.  From  a  comparison  of  C=€  bond  lengths 
with  those  of  ethylene  (1.333  A)  and  nitroethylene  (1.329  A),  it  appears  that  the  nitro  group 
strengthens  the  olefinic  bond  whereas  it  is  weakened  by  the  additional  substituents  (Hi,  NH„ 
and  to  a  smaller  degree  CH,.  This  effect  may  reflect  extended  conjugation. 
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The  energy  difference  between  the  2-nitrovinyl  alcohol  and  aci-nitroacetaldehyde 
tautomers  is  only  4.8  kcal/mol  with  a  small  barrier  of  S.8  kcal/mol  for  the  forward  reaction. 

A  significantly  larger  energy  difference  of  14.5  kcal/mol  is  calculated  for  the  2-nitrovinyl 
amine  and  oci-nitroacetaldimine  tautomers,  of  which  the  latter  has  a  kinetic  stability  of  only 
2.2  kcal/mol  at  MP4/6-31G*  but  none  when  zero-point  energy  corrections  are  included.  1- 
Nitropropylene  is  18.7  kcal/mol  more  stable  than  its  nitronic  acid  and  would  minimally 
require  a  large  42.4  kcal/mol  for  such  a  rearrangement  to  occur.  Figure  9  illustrates  not  only 
the  large  differences  between  these  isoelectronic  systems  but  also  lists  the  significant 
dependency  of  energy  differences  on  the  theoretical  levels  employed.  The  tendency  for 
tautomerism  appears  to  decrease  both  kinetically  and  thermodynamically  on  introducing  the 
substituents  in  the  order  OH,  NH^  and  CH3.  It  is  recognized  that  the  MP4/6-31G*  level  of 
theory  is  not  adequate  for  a  final  evaluation  of  the  potential  energy  surface  (see  Section  II) 
and  hence  these  energy  differences  are  likely  to  reduce  at  more  sophisticated  levels. 

Competing  tautomers  are  those  resulting  from  formal  [1,3]-H  shifts  and  these  are  2- 
nitroacetaldehyde,  2-nitroacetaldimine,  and  3-nitropropene.  Their  MP4/6-31G*  (+  ZPE) 
energy  differences  with  the  nitroethylene  tautomers  are  -  3.6,  +2.4,  and  +1.1  kcal/mol 
respectively.  This  means  that  in  these  systems  the  keto  enol  energy  differences  amounts  to 
only  4.0  kcal/mol,  whereas  it  is  17.2  kcal/mol  for  acetaldehyde  vinyl  alcohol  at  the  same 
level  of  theory.  Considering  that  this  latter  value  reduces  to  10.9  kcal/mol  at  Gl,  and 
applying  this  'correction*  also  to  the  conjugated  analogue  leads  to  the  likely  conclusion  that  2- 
nitrovinyl  alcohol  is  the  global  minimum.  The  same  argument  applies  to  the  other  compounds 
and  these  are  summarized  in  Figure  10. 
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The  Effects  of  Electron  Correlation  on  Geometries 
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Figure  5.  Structures  of  the  Tautomers  of  Nitrovinyl  Alcohol. 


Figure  7.  MP2/6-31G*  Structures  of  the  cis  Tautomers  of  2-Nitrovinyl  Amine  and  their 
Transition  Structures. 
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Structure. 
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TABLE  8.  ABSOLUTE  ENERGIES  (in  -  hartrees)  OF  SUBSTITUTED  NITROETHYLENE 

TAUTOMERS  AND  THEIR  ANIONS.* 


Compound 

HF/3-21G 

HF/6-31G* 

MP2/6-31G* 

MP4/6-31G** 

C^NO, 

16  C, 

354.34455  (0) 

356.37129  (0) 

357.36069  (0) 

357398291 

16a 

C. 

354.31813  (0) 

356.35286  (0) 

357.34116  (0) 

35734189 

17 

c. 

354.33883  (0) 

356.35703  (0) 

357.35329  (0) 

357.38967 

18 

c. 

354.33310  (1) 

356.34804  (1) 

357.34911  (1) 

357.38392 

19 

C, 

354.33225  (0) 

356.37765  (0) 

357.36617  (0) 

357.40566 

TO 

c. 

353.78110  (0) 

355.81821  (0) 

356.80925  (0) 

356.84379 

CW*A 

21  C. 

334.63553  (0) 

336.54755  (0) 

337.52181  (0) 

337.56267 

22 

c. 

334.61206  (0) 

336.51574  (0) 

337.49950  (0) 

337.53963 

23 

c. 

334.61017  (1) 

336.50763  (1) 

337.49748  (1) 

337.53609 

24 

c, 

334.60560  (0) 

336.54003  (0) 

337.51459  (0) 

337.55875 

25 

c. 

334.02769  (0) 

335.95711  (0) 

336.93604  (0) 

336.97493 

CftNOj 

26  C. 

318.71225  (0) 

320.54169  (0) 

321.48544  (0) 

32133231 

26a 

c. 

318.71184  (0) 

320.54120  (1) 

321.48438  (1) 

321.53118 

27 

c, 

318.71557  (0) 

320.51239  (0) 

321.45489  (0) 

32130247 

28 

c» 

318.64118  (1) 

320.44925  (1) 

321.41821  (1) 

321.45982 

29 

c, 

318.70877  (0) 

320.54121  (0) 

321.48403  (0) 

321.53107 

30 

c. 

318.13774  (0) 

319.95544  (0) 

320.90327  (0) 

320.94546 

*  Values  in  parentheses  are  the  number  of  imaginary  frequencies. b  MP2/6-31G* 
geometries. 


23 


"  BLE  9.  RELATIVE  ENERGIES  (in  kcal/inol)  OF  SUBSTITUTED  NITROETHYLENE 

TAUTOMERS  AND  THEIR  ANIONS. 


Compound 

HF/3-21G 

HF/6-31G* 

MP2/6-31G* 

MP4/6-31G** 

MP4/6-31G** 

+  ZPE 

CjHsNO, 

16 

c. 

0.0 

0.0 

0.0 

0.0 

0.0 

16a 

c. 

16.6 

11.6 

12.3 

12.0 

12.4 

17 

c. 

3.6 

8.9 

4.7 

5.4 

4.8 

18 

c. 

7.2 

14.0 

7  3 

9.0 

5.8 

19 

Ci 

7.7 

-4.0 

-3.4 

-4.6 

-3.6 

20 

C, 

353.6 

347.1 

346.0 

348.0 

339.1 

(W*A 

21 

C. 

0.0 

0.0 

0.0 

0.0 

0.0 

22 

Q 

14.7 

19.9 

14.0 

14.5 

14.5 

23 

c. 

15.9 

25.1 

15.3 

16.7 

13.9 

24 

ct 

18.8 

4.7 

4.5 

2.5 

2.4 

25 

C. 

381.4 

370.5 

367.6 

368.8 

360.4 

QHjNO* 

26 

C. 

0.0 

0.0 

0.0 

0.0 

0.0 

26a 

C. 

0.3 

0.3 

0.7 

0.7 

0.6 

27 

C, 

-2.1 

18.4 

19.2 

18.7 

18.7 

28 

C, 

44.6 

58.0 

42.2 

45.5 

42.4 

29 

C, 

2.2 

0.3 

0.9 

0.8 

1.1 

30 

C. 

360.5 

367.8 

365.3 

368.3 

360.0 

*  MP2/6-31G*  geometries. 
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Figure  10.  MP4/6-31G*  Energy  Comparisons  On  kcal/moi)  between  Substituted  Nitroethylenes 
Tautomers  and  the  Parent  Nitro/fld-Nitro,  Keto/Enol,  and  Imine/Enamine  Tautomers. 
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Bond  Critical  Point  Data  for  Substituted  Nitroethylenes. 
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Figure  11.  Hydrogen  Bonding  at  MP2/6-31G*  for  Substituted  Nitroethyienes  and  N-OH  Bond  Strengths  in 
their  aci-Nitro  Tautomers. 


SECTION  IV 


TAUTOMERISM  IN  SUBSTITUTED  NITROBENZENES 

1.  INTRODUCTION 

We  consider  o-nitrophenol,  o-nitroaniline,  ami  o-nitrotoluene  in  this  section.  Their 
HF/3-21G  and  HF/6-31G*  optimized  structures  and  those  of  their  ori-tautomers  and  the 
transition  structures  for  1,5-H  shifts  are  displayed  in  Figures  13-15.  Total  HF  and  MP/6- 
31G*  and  relative  energies  are  given  in  Table  10,  and  visualized  in  Figures  16  and  17. 

2.  BACKGROUND 

The  study  by  Politzer  and  cowodcets  on  the  tautomerism  of  o-nitrophenoi  (Reference 
1),  and  their  data  on  this  tautomeric  pair,  and  their  unsuccessful  attempts  of  finding 
tautomerism  in  other  substituted  nitroaromatic  compounds*4  form  the  underpinnings  of  this 
investigation.  Of  particular  interest  are  the  reported  tautomer  of  o-nitrophenol  and  the  author's 
inability  to  locate  similar  tautomeric  structures  for  o-nitroaniline  and  o-nitrotoluene.  These 
calculations  were  performed  at  die  3-2 1G  level  of  theory  and  the  results  are  counterintuitive 
to  those  in  die  previous  section.  It  is  important  to  stress  that  the  conformation  of  aci-2- 
nitrophenol,  reported  by  die  Politzer  group,  has  a  syn  orientation  of  die  (O)NO — H  bond  with 
respect  to  die  neighboring  C=0  group  with  an  energy  difference  of  15.4  kcaMmol  between  the 
two  tautomers. 

3.  DISCUSSION 

The  HF/6-31G*  optimized  structures  of  the  nitro  compounds  are  in  reasonable 
agreement  with  X-ray  structures  (Reference  42).  The  largest  differences  are  found  for  die 
C-N  and  N-O  bonds  as  already  discussed  in  Section  IL  The  (0)N0.~H0  distance  in  o- 
nitrophenol  is  relatively  short  with  1.826  A  at  HF/6-31G*  (1.91  A,  X-ray;  1.748  A,  HF/3- 
21G)  and  this  is  also  the  case  for  die  (0)N0.~HNH  distance  in  o-nitroaniline  for  which  the 
HF/6-31G*  separation  is  calculated  at  1.991  A  (1.955  A,  X-ray).  However,  numerous 
attempts  to  calculate  die  nitronic  add  tautomer  of  both  o-nitrophenol,  (as  reported  by  Politzer 
and  coworkers),  and  o-nitroaniline  at  HF/3-21G  proved  unsuccessful  These  tautomers  could 
only  be  calculated  with  the  larger  6-3 1G*  basis  set  at  which  level  also  die  transition  structures 
were  determined  and  characterized.  However,  die  kinetic  stability  of  both  nitronic  acids  is 
less  than  1  kcal/mol  Single  point  calculations  that  indude  the  effects  of  electron  correlation 
at  MP2/6-31G*  give  lower  relative  energies  for  the  transition  structures  than  for  the  nitronic 
add  tautomers.  This  is  illustrated  in  Figure  17  and  strongly  suggests  that  the  nitronic  add 
tautomers  of  o-nitrophenol  and  o-nitroaniline  are  not  likely  to  have  a  trans  oci-nitro 
conformation.  In  this  context  it  is  relevant  to  recall  that  die  oci-nitro  tautomer  of  1- 
nitropropene  has  a  cis  conformation  of  die  (O)NO — H  group. 
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Optimization  of  the  nitronic  adds  of  o-nitrophenol,  o-nitroaniline,  and  o-nitrotoluene 
with  a  ds  orientation  within  the  0«-N — O — H  oci-nitro  group  results  in  local  minima  in  all 
cases.  The  transition  structures  for  the  [1,5]-H  shift  were  also  identified  and  the  results  are 
displayed  in  Figure  17.  Three  observations  can  be  made.  Firstly,  the  nitronic  adds  have 
significant  kinetic  stabilities,  Le.  14.0  kcal/mol  for  o-nitrophenol,  4.2  kcal/mol  for  o-nitro- 
anfline,  and  6.2  kcal/mol  for  o-nitrotoluene  (MP2/6-31G*  +  ZPE).  Secondly,  the  energy 
differences  (same  level)  between  the  tautomers  increase  from  o-nitrophenol  (30.7  kcal/mol)  to 
o-nitroaniline  (43.3  kcal/mol)  to  o-nitrotoluene  (46.3  kcal/mol).  Thirdly,  these  kinetic  and 
thermodynamic  energy  differences  are  much  larger  than  those  of  the  substituted  nitroethylenes 
(except  for  nitropropylene)  and  those  of  the  O-N — O — H  trans  conformations  of  die 
aromatic  nitronic  adds.  The  influence  of  electron  correlation  (at  MP2)  on  the  energies  is  only 
modest  except  for  die  rearrangement  barrier  of  o-nitrotoluene.  It  must  be  emphasized  dial 
whereas  all  nitro  compounds  show  small  distortions  from  planarity  at  HF/6-31G*  the 
transition  structures  deviate  significantly  from  planarity,  both  for  the  transferring  H  and  for 
die  aromatic  ring.  This  is  illustrated  in  Figures  13-16.  The  main  conclusion  from  this  section 
is  that  nitro  «*  aci- nitro  tantomerism  is  feasible  for  o-substituted  nitrobenzenes  only  when  the 
nitronic  add  has  a  ds  CH-N — O — H  conformation. 


30 


NHro  **  acbNltro  Tautomerism  in  ortho-Nitrophenoi 
Nftro  TS  acLNitro 


AE  (MP2/6-31G*):  45.3  kcat/mol  32.4  kcal/mol 


hf/wiq* 

(HF/3-21G) 


Figure  13.  Tautomers  of  o-Nitrophenol. 
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Nltro  <■*  acl-Nltro  Tautomerlsm  In  ortho-NItrotolueno 


Figure  15.  Tautomers  of  o-Nitrotoluene. 
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Energies  of  Substituted  Nltroaromatlc  Tautomers 

(tram  icWtfO,  O-N-O-H) 

MP2/6-31G* 
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Figure  16.  Energy  Profile  of  Tautomerism  of  o-Nitrophenil  and  o-Nitroaniline  to  trans  Nitronic  Acids. 
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SECTION  V 


TAUTOMERISM  IN  SUBSTITUTED  DINTTROBENZENES 

1.  INTRODUCTION 

We  consider  2,6-dinitrophenol,  2,6-dinitroaniline,  and  2,6-dinitrotoluene.  The  HP/3- 
210  and  HF/6-31G*  optimized  parameters  of  the  nitro,  aci- nitro,  and  related  transition 
structures  are  shown  in  Figures  18-20.  Total  and  relative  3-21G  and  6-3 1G*  energies  are 
given  in  Table  11 

2.  DISCUSSION 

Many  of  the  same  observations  as  those  discussed  in  the  section  on  substituted 
nitrobenzenes  also  apply  to  die  dinitrobenzenes.  However,  the  buttressing  effect  of  die  second 
nitro  group  is  evident  ami  this  affects  the  three  nitro  compounds  in  different  ways.  One  nitro 
group  is  significantly  rotated  (37*)  in  2,6-dinitrophenol,  both  nitro  groups  are  planar  in  die 
aniline  derivative  due  to  H-bonding  ((0)N0..JiNH  =  1.959  A),  while  both  nitro  groups  are 
rotated  (35*)  in  the  toluene  derivative  to  minimize  the  repulsion  with  the  methyl  group.  The 
same  applies  to  the  nitronic  acids  with  NOj  bond  rotations  of  35*  for  the  dinitrophenol 
tautomer  and  23*  for  the  dinitrotoluene  tautomer,  and  die  dinitroaniline  tautomer  is  now  also 
non-planar.  All  tautomers  are  minima.  The  transition  structures  for  the  [1,5]-H  transfer  were 
also  characterized  at  both  HF/3-21G  and  HF/6-31G*.  Comparison  of  the  HF/6-31G*  energy 
differences  (A£)  between  the  nitro  -  oci-nitro  tautomeric  pairs,  (barriers  are  in  parentheses), 
with  those  of  the  mononitro  species  highlight  subtle  but  important  differences.  For  example, 

A E  for  2,6-dinitrophenol  is  29  (34)  kcal/mol  whereas  a  much  larger  barrier  of  46  kcal/mol 
(A E  -  31)  is  calculated  for  o-nitrophenoL  The  2,6-dinitro  derivative  of  aniline  has  a  much 
larger  A E  of  45  (49)  kcal/mol  than  that  of  phenol,  but  in  contrast  it  compares  with  the  43  (48) 
kcal/mol  calculated  for  o-nitroaniline.  The  aniline  nitro  oc/'-nitro  energy  differences,  but  not 
the  barriers,  are  comparable  in  magnitude  to  the  41  (65)  kcal/mol  for  2,6-dinitrotoluene  and 
the  43  (69)  kcal/mol  for  o-nitrotoluene.  The  effects  of  electron  correlation  in  the  AE  values 
could  not  be  evaluated  directly  due  to  limitations  in  resources.  However,  if  the  influence  of 
electron  correlation  (MP2)  on  AE  for  the  dinitro  derivatives  is  similar  to  that  for  the 
mononitro  compounds,  then  it  must  be  assumed  that  the  barrier  for  [1,5]-H  transfer  will 
reduce  significandy  at  correlated  levels.  Nevertheless,  it  appears  that  the  tautomeric  AE  is 
slightly  less  for  dinitrotoluene  than  it  is  for  dinitroaniline. 
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Tautomers  of  2,6-DinHrophonoi 


Tautomers  of  2,6-DinHroanMne 


HF/6-31G* 

(HF73-21G) 


AE:  46.7  kcal/mol 
(41.7) 


barrier  49.4  kcal/mol 
(46.9) 


1.197 

OJ*l 


Figure  19.  Optimized  Tautomeric  Structures  of  2,6-Dinitroanilme. 
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Toulomon  of  2,6-UnHrotokmto 
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AE:  41.3  kcal/mol 
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Figure  20.  Optimized  Tautomeric  Structures  of  2,6-Dinitrotohieiie. 


SECTION  VI 


TAUTOMERKM  IN  PICRIC  ACID,  PICRAMIDE,  AND  TNT 

1.  INTRODUCTION 

la  this  lection  we  consider  picric  add,  picramide,  and  TNT.  Their  folly  optimised  and 
characterized,  noo-symmetric  duo,  ori-nitro,  and  corresponding  transition  structures  are 
displayed  hi  Figures  21-23.  Total  and  relative  HEV3-21G  and  HF/6-31G*  energies  are  given 
inlhhle  12. 

2.  DISCUSSION 

The  optimized  nitro  compound  structures  compare  well  to  those  reported  from  X-ray 
structure  determinations  (Reference  43).  As  referred  to  in  the  previous  sections,  the  largest 
discrepancies  are  with  the  N-O  and  C-N  bonds  but  in  these  cases  also  deviations  are  found 
for  the  rotation  of  the  nitro  groups.  As  compared  to  the  dinitro  derivatives  of  the  previous 
section,  addition  of  the  third  nitro  group  has  s  subtle  but  important  influence  beyond  that  of  a 
spectator  at  the  4-positkm.  The  geometrical  parameters  for  picric  add  and  2,6-dinitrophcnol 
are  very  similar,  even  for  the  o-nitro  group  HW6-31G*  bond-rotations  of  35*  and  37\ 
respectively,  wife  corresponding  rotations  of  38*  and  33*  for  their  nitronic  adds.  However, 
tiie  two  X-ray  crystal  stracturea  of  picric  add  indicate  a  much  smaller  rotation  of  '0111/  17- 
20*,  which  interestingly  compares  well  with  the  HF/3-21G  value  of  20*.  Picramide  has  a 
planar  structure  which  compares  well  the  X-ray  structure  and  is  similar  to  that  of  2jth 
dinitroaniline.  This  suggests  that  it  is  also  similarly  stabilized  by  two  H-bonds  between  the 
nitro  groups  and  the  amine  group.  TNT  is  noo-planar  but  does  have  a  plane  of  symmetry. 

I  Ac  2,6-dimtrotohiene,  its  nitro  groups  are  rotated  (36*,  6-31G*;  43-31,  exp.)  with  one  of  the 
oxygens  of  each  nitro  group  bisecting  the  neighboring  methyl  group.  This  conformation  is  the 
same  as  for  the  mono-  and  dinitro  derivatives.  The  nitronic  acid  tautomer  of  TNT  has  s 
'benzene*  ring  that  is  significantly  distorted  from  planarity.  This  is  likely  the  result  of  steric 
congestion  between  the  methylene  group  and  its  neighboring  oci-nitro  and  nitro  groups.  The 
nitro  group  rotation  is  23*. 

The  energy  differences  between  the  nitro  compounds  and  their  oci-nitro  tautomers 
amplify  the  trend  already  observed  for  die  dinitro  compounds.  Barriers  for  [1,5]-H  transfer 
from  the  nitro  to  the  nitronic  acid  tautomers  are  given  in  parentheses.  Thus,  A E  (HF/6-31G*) 
for  picric  add  has  reduced  to  23  (34)  from  29  (34)  kcal/mol  for  dinitrophenoL  likewise,  A E 
for  TNT  of  39  (63)  kcal/mol  is  smaller  than  the  41  (63)  kcal/mol  for  dinitrotoloene.  In 
contrast,  a  larger  A£  of  48  (31)  kcal/mol  is  calculated  for  picramide  as  compared  to  a  AE  of 
43  (49)  for  dinitroaniline.  Even  though  these  changes  are  small  and  can  be  individually 
explained  by  shortcomings  in  the  theoretical  levels  employed,  they  continue  the  trend  already 
observed  for  die  dinitro  derivatives. 
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3.  CORRELATION  BETWEEN  SENSITIVITIES  AND  TAUTOMERISM 

Important  links  between  impact  and  shock  sensitivities  and  the  accessibility  of 
tautomerism  may  be  explored  even  though  the  levels  of  theory  used  for  these  large  systems 
are  inadequate  to  make  definitive  conclusions.  Firstly,  the  energy  difference  between  the 
tautomers  of  picric  arid  (25  kcal/mol,  HF/6-31G*)  as  well  as  the  barrier  linking  them  (34 
kcal/mol)  is  by  far  the  smallest  of  the  set  of  nitroaromatics  studied.  Picric  arid  is  also  die 
most  sensitive  compound.  The  calculated  relative  energy  differences  of  both  the  nitronic  arid 
and  the  transition  structure  are  expected  to  be  reduced  further  by  the  effects  of  electron 
correlation  beyond  MP2.  Secondly,  die  energy  difference  between  the  tautomers  of  picramide 
(48  kcal/mol,  HF/6-31G*)  is  somewhat  larger  than  that  between  die  TNT  tautomers  (39 
kcal/mol).  Picramide  is  also  a  less  sensitive  compound  than  TNT.  Both  picramide  and  TNT 
have  significant  barriers  for  tautomerism  (51  and  65  kcal/mol,  respectively,  HF/6-31G*),  but 
these  are  expected  to  reduce  somewhat  when  the  effects  of  electron  correlation  are  taken  into 
account  with  the  larger  effect  expected  for  TNT.  Thirdly,  when  the  MP2/6-31G*  reduction  in 
the  barrier  for  [1,5]-H  transfer  for  o-nitrotoluene  is  applied  to  TNT,  an  activation  barrier  of  52 
kcal/mol  is  estimated.  This  barrier  is  very  similar  to  the  experimental  induction  activation 
energy  for  thermal  decomposition  of  liquid  TNT.  Fourthly,  the  tautomerism  barriers  for 
picramide  and  TNT  are  significantly  higher  than  that  of  picric  arid  and  may  suggest  alterna¬ 
tive  decomposition  pathways  for  which  there  is  experimental  support 
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Figure  21.  Optimized  Tautomeric  Structures  of  Picric  Acid. 
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Figure  22.  Optimized  Tautomeric  Structures  of  Picricamide. 
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Figure  23.  Optimized  Tautomeric  Structures  of  TNT 


TABLE  12.  TOTAL  (in  -  hartrees)  AND  RELATIVE  (in  kcal/mol)  ENERGIES  OF  SUB! 
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SECTION  vn 


CYCUZATIONS  REACTIONS  OF  SUBSTITUTED  NITROETHYLENES 

1.  INTRODUCTION 

In  the  context  of  unraveling  the  intricacies  that  underlie  the  combustion  mechanisms  of 
less  sensitive  nitroaromatic  explosives,  it  is  of  importance  to  determine  the  key  intermediates 
that  would  trigger  such  a  process.  The  charge  of  this  contract  research  was  to  establish 
whether  or  not  nitronic  acids  are  viable  candidates  to  occupy  such  a  central  role.  This  was 
answered  in  the  above  sections.  Here  logical  next  steps  are  explored.  If  the  nitronic  acids 
play  indeed  a  pivotal  role  it  follows  that  their  decomposition  pathways  must  be  energetically 
competitive  with  those  that  involve  the  nitroaromatics  directly.  The  two  processes  that  were 
further  investigated  are  bond  dissociation  energies  and  cyclization  reactions  to  alternative 
products.  These  studies  were  conducted  beyond  the  contract  charge  to  probe  future  directions 
to  determine  in  greater  detail  the  combustion  mechanism  of  polynitroaromatics. 

2.  BOND  DISSOCIATION  ENERGIES 

A  large  spectrum  of  radicals  that  constitute  fragments  of  simple  nitroalkanes  and 
substituted  nitroalkenes  was  calculated.  These  were  all  fully  optimized  and  characterized  at 
various  levels  of  theory.  Table  13  summarizes  the  energy  data  and  also  lists  the  <S*,1>  data 
which  underscore  that  not  all  radicals  were  calculated  satisfactorily  as  identified  by  their 
significant  deviations  from  0.7S.  It  is  recognized  that  more  sophisticated  levels  of  theory  (e.g. 
MRCI)  are  required  to  accurately  calculate  bond  dissociation  energies  in  particular  for  the 
nitroso  radicals  for  which  various  electronic  states  can  be  envisioned.  Nevertheless,  some 
insight  may  be  obtained  from  the  listed  data.  Several  combinations  can  be  made  to  obtain 
crude  estimates  for  various  bond  dissociation  energies.  We  concentrate  only  on  the  nitro  and 
nitronic  acid  groups  as  summarized  in  Tables  14-17.  The  C-N02  bond  dissociation  energies 
of  nitromethane  and  nitroe thane  are  remarkably  similar  to  the  experimentally  determined 
values.  The  C-N02  bond  dissociation  energy  increases  by  introducing  a  double  bond  that  is 
conjugated  to  the  nitro  group,  i.e.  nitroethylene.  Even  larger  values  are  obtained  for  the 
hydroxy  and  amino  substituted  olefins  (see  Table  14).  These  data  suggest  that  in  the  absence 
of  steric  effects  the  C-N02  bond  is  strengthened  by  conjugative  effects. 

The  (O)N-OH  bond  dissociation  energies  of  the  nitronic  acids  of  the  substituted 
nitroethylenes,  listed  in  Table  17,  underscore  that  inadequate  levels  of  theory  have  been 
employed.  The  data  do  suggest  that  PMP2/6-31G*  overestimates  the  bonding  energies  and 
that  the  N-OH  bond  is  likely  to  be  weak,  but  this  must  be  verified  at  more  sophisticated 
levels  of  theory.  It  is  plausible,  however,  that  the  nitronic  acids  are  likely  candidates  to 
generate  nitroso  compounds,  which  are  experimentally  abundant  products  in  the 
decompositions  of  nitroaromatics. 
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TABLE  13.  ABSOLUTE  ENERGIES  (is  hmw)  OP  NITRO METHANE,  NITRO  ETHANE, 
NITRO  ETHYLENE,  SUBSTITUTED  NITRO  ETHYLENE,  AND  THEIR  FRAGMENTS. 


Molecule 

Sym.  UHF/3-21G  UHF/6-31G* 

PMP2/6-31G* 

PMP4/6-31G** 

ZPE*  I* 

s1' 

CH3NO3 

c. 

-242.25586 

-243.66199 

-244.34534 

-24437427 

30.80 

0 

■ 

C3H3NO2 

c. 

-279.88906 

-281.30409 

-282.31131 

-282.34769 

34.00 

0 

- 

CjHjNOj 

c. 

-281.08265 

-282.70109 

-283.52047 

-283.55884 

48.28 

0 

- 

HO-CH=CH-NOj 

c. 

-354.34455 

-356.37129 

-35736070 

-357.39829 

37.85 

0 

- 

H3N-CH=CH-NOj  c. 

-334.63553 

-336.54755 

-33732181 

-337.56267 

44.86 

0 

- 

HsC-CH=CH-N03  C, 

-318.71225 

-320.54169 

-321.48544 

-321.53231 

51.15 

0 

- 

0=CH-CH=N0aH  C, 

-354.33883 

-356.35703 

-35735329 

-357.38967 

37.29 

0 

- 

HN=CH-CH=N03H  c. 

-334.61206 

-336.51574 

-337.49950 

-337.53963 

44.87 

0 

- 

h2c=ch-ch=no2h  c. 

-318.71557 

-320.51239 

-321.45489 

-321.50247 

51.12 

0 

- 

HO* 

c. 

-74.97023 

-75.38228 

-75.52475 

-7533695 

5.15 

0 

0.75 

NOj* 

C, 

-202.8435 

-204.03149 

-204.57253 

-20439027 

5.56 

0 

0.76 

ch3* 

c. 

-39.34261 

-39.55899 

-39.67504 

-39.69056 

17.50 

0 

0.75 

CHO* 

c. 

-112.60380 

-113.47660 

-113.54277 

-113.55751 

8.12 

0 

0.76 

CHjN* 

c. 

-92.86450 

-93.39296 

-93.66420 

-93.68376 

14.67 

0 

0.79 

CHjNOa* 

C. 

-241.61810 

-242.98587 

-243.62084 

-243.64652 

21.53 

0 

1.12 

C3HS* 

Cs 

-76.96280 

-77.39029 

-77.62332 

-77.64605 

21.85 

0 

0.89 

c3h5* 

c. 

-78.16365 

-78.59715 

-78.84679 

-78.87368 

35.70 

0 

0.75 

C3H30: 

c. 

-150.72351 

-151.58631 

-151.99331 

-152.02453 

17.71 

0 

- 

C3H,N: 

C, 

-131.00487 

-131.75694 

-132.15135 

-132.18727 

25.54 

0 

- 

c3hso* 

c. 

-151.39772 

-152.24291 

-152.65714 

-152.68242 

25.11 

0 

0.87 

c3h4n* 

c. 

-131.68236 

-132.41496 

-132.81670 

-132.84628 

32.44 

0 

0.88 

c3h3no3* 

c. 

-279.38148 

-280.96491 

-281.72012 

-281.75418 

24.41 

0 

0.99 

c3h3n3o* 

Ci 

-259.63907 

-261.11555 

-261.85738 

-261.89768 

33.07 

0 

0.89 

C3H4: 

c. 

-115.12652 

-115.77656 

-116.14060 

-116.17728 

32.62 

0 

- 

c3h5* 

c. 

-115.78488 

-116.42915 

-116.79742 

-116.83097 

39.47 

0 

0.88 

CsH4NO* 

c. 

-243.76104 

-245.13364 

-245.83927 

-245.88786 

38.03 

0 

1.15 

H03N:(trans) 

c. 

-203.46601 

-204.63768 

-205.17335 

-205.19587 

13.03 

0 

- 

H03N:(cis) 

c. 

-203.46850 

-204.63994 

-205.17536 

-205.19787 

13.08 

0 

- 

a  PMP4(SDTQ)/6-31G*//MP2/6-31G*.  *  Zero-point  energies  scaled  by  0.9. 
*  Number  of  imaginary  frequencies.  d  After  annihilation  at  PMP2/6-31G* 
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TABLE  14.  C-N02  BOND  DISSOCIATION  ENERGIES  (in  kcal/mol)  FOR 
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Relative  energies  in  the  second  column  are  corrected  for  zero-point  energies. 


TABLE  16.  C=N  BOND  DISSOCIATION  ENERGIES  (in  kcal/moi)  FOR  SUBSTITUTED  ari-NITROETHYL 
Compound  UHF/3-21G  UHF/6-31G*  PMP2/6-31G*  PMP4/6-31G** 
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Relative  energies  in  the  second  column  are  corrected  for  zero-point  energies. 


3. 


CYCUZATTON  REACTIONS 


The  objective  of  the  model  study  described  in  this  section  was  to  elaborate  a  pathway 
that  could  mimic  possible  formation  of  anthranilic  acid  from  nitroaromatics  by  means  of 
nitronic  acids.  Such  pathways  can  be  elaborated  by  both  radical  and  acid/base  mechanisms. 

In  either  case,  the  involvement  of  nitronic  acids  would  be  an  important  component 
Therefore,  valuable  information  with  respect  to  the  integrity  of  nitronic  acids  was  sought. 

We  considered  the  lowest  energy  pathway  for  cyclization  of  nitro vinyl  alcohol  (16), 
nitrovinyl  amine  (21),  and  3-nitropropene  (26)  via  non-ionic  and  non-radical  routes.  The  total 
and  relative  HF/3-21G,  HF/6-31G*,  and  MP2/6-31G*  energies  of  the  various  feasible 
intermediates  and  associated  transition  structures  are  listed  in  Table  18-20,  respectively.  The 
pathways  are  depicted  schematically  in  the  corresponding  Figures  24-26,  which  display  the 
structures  and  the  MP2  energy  profiles  for  the  cyclization  reactions.  There  is  some  necessary 
overlap  with  Section  IIL 

Cyclization  is  feasible  for  2-nitro vinyl  amine  and  1-nitropropene  whose  products  can 
undergo  dehydration,  which  in  both  cases  renders  the  entire  reaction  scheme  exothermic  at  the 
correlated  level  of  theory.  This  exothermicity  is  largest  for  1-nitropropene,  which  also  has  the 
lower  energy  barriers  for  the  primary  cyclization  reaction,  Le.  49.7  kcal/mol  vs  66.3  kcal/mol 
relative  to  21  and  26,  respectively.  No  cyclization  product  could  be  found  at  the  correlated 
level  for  2-nitrovinyl  alcohol,  (at  HF/6-31G*  a  cyclic  structure  was  obtained  with  an 
associated  barrier  for  formation  of  88.3  kcal/mol,  relative  to  16).  While  it  is  recognized  that 
the  peroxy  group  is  theoretically  difficult  to  calculate  accurately,  the  large  endothermicity 
makes  cyclization  an  unlikely  prospect  This  is  also  the  case  for  the  cyclization  of  2- 
nitro vinyl  amine  which  is  endothermic  by  61.3  kcal/mol.  However,  the  cyclization  product 
and  the  nitronic  adds  of  1-nitropropene  are  nearly  of  the  same  energy  and  only  13-20 
kcal/mol  less  stable  than  26,  with  corresponding  similar  transformation  barriers  of  42-30 
kcal/mol. 

The  reaction  profiles  istrate  that  cyclization  of  1-nitropropene  is  feasible  and  this 
compares  well  with  the  experimentally  observed  anthranilic  adds.  Cyclization  of  2-nitrovinyl 
amine  is  possible  but  may  be  energetically  prohibitive  wi'h  little  kinetic  stability  of  the 
product  Cyclization  of  2-nitrovinyl  alcohol  appears  not  to  be  feasible. 

Because  of  the  relation  with  the  nitroaromatics  it  is  of  interest  to  briefly  compare  the 
various  nitronic  add  conformations.  As  demonstrated  in  Sections  IV-VI  all  aromatic  nitronic 
adds  have  a  cis  conformation  of  the  Ck-N-O-H  group;  it  was  shown  in  Section  IV,  Table  10 
and  Figures  13  and  16,  that  the  trans  conformations  have  no  kinetic  stability.  In  the  case  of 
the  ethylenes  the  kinetic  stability  of  the  cis  forms  (32  and  38)  are  8-10  kcal/mol  less  stable 
(barrier  ca.  13  kcal/mol)  than  the  corresponding  trans  forms  (17  and  22).  Interestingly,  in  all 
three  systems  (2-nitrovinyl  alcohol,  2-nitrovinyl  amine,  and  1-nitropropene)  the  nitronic  adds 
with  the  N-OH  farthest  from  the  substituent  (34,  40,  and  45)  are  the  most  stable  cis 
0«-N-0-H  conformers.  These  have  also  significant  kinetic  stability,  although  under 
experimental  conditions  (liquid,  solid)  such  transformations  may  occur  ionically  (Section  II). 

It  is  then  relevant  to  note  that  these  nitronic  adds  34, 40,  and  45  are  'only'  11.3,  22.4,  and 
1 3.0  kcal/mol  less  stable  than  their  corresponding  most  stable  tautomers.  Especially  the 
relative  order  is  of  interest  as  it  compares  well  with  that  calculated  for  picric  add,  picramide, 
and  TNT. 
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TABLE  16.  TOTAL  (in  hnitnan)  AMD  RELATIVE  (is  kcnl/anl)  ENERGIES  OF  INTERMEDIATES 
AND  TRANSITION  STRUCTURES  FOR  THE  CYCUZATION  OF  3-N1TROETHENOL. 


Structures 

Sym. 

HF/3-21G 

mbs. 

ral. 

HF/6-31G 

abs. 

» 

rel. 

MP2/6-31G 

aba. 

* 

rel. 

C3H3NO3 

16 

c. 

—354.34455 

0.0 

-356.37129 

0.0 

-357.36070 

03 

18  (16  -♦  17) 

c. 

—354.33310 

7.2 

-35634804 

14.0 

-357.34911 

73 

17 

c. 

-35433883 

3.6 

-35635703 

8.9 

-357.35329 

4.7 

SI  (17  -»  32) 

c. 

-35432119 

14.2 

-35633876 

20.4 

-35733133 

18.4 

32 

c. 

-35433333 

7.0 

-356.34941 

13.7 

-357.34040 

12.7 

33  (32  -»  34) 

c. 

-35439237 

32.7 

-356.29333 

48.9 

-357.29547 

40.9 

34 

C. 

-35433411 

6.5 

-35634947 

13.7 

-357.34262 

113 

35  (34  -►  36) 

Cl 

-35434547 

62.2 

-356.23025 

88.5 

36 

Cl 

-35439131 

33.4 

-35635689 

713 

TABLE  19.  TOTAL  (in  knitnaa)  AND  RELATIVE  (in  kcnl/nwl)  ENERGIES  OF  INTERMEDIATLS 
AND  TRANSITION  STRUCTURES  FOR  THE  CYCLIZATION  OF  S-NTTROVINYLAMINE. 


Structures 

Sym. 

HF/3-21G 

abs. 

rd. 

HF/6-31G* 

abs. 

rel. 

MP2/6-31G 

abs. 

* 

rel. 

C3H4N3O3 

21 

c. 

-334.63553 

0.0 

-336.54755 

0.0 

-337.52181 

0.0 

23  (21  ->  22) 

c. 

-334.61017 

15.9 

-336.50763 

25.1 

-337.49748 

15.3 

22 

c. 

-334.61206 

14.7 

-336.51574 

19.9 

-337.49950 

14.0 

37  (22  —  38) 

c. 

-334.59238 

27.1 

-336.49782 

31.2 

-337.47537 

29.1 

38 

c. 

-334.60233 

20.8 

-336.50685 

23.5 

-337.48295 

24.4 

39  (38  —  40) 

c. 

-334.56441 

44.6 

-336.45321 

59.2 

-337.43985 

51.4 

40 

c. 

-334.60294 

20.5 

-336.50700 

25.5 

-337.48612 

22.4 

41  (40  -+  42) 

Ci 

-334.52723 

67.9 

-336.40669 

88.4 

-337.41616 

66.3 

42 

Ci 

-334.58858 

29.5 

-336.45978 

55.1 

-337.42408 

61.3 

43 

c. 

-259.06708 

-260.54385 

-261.33605 

H3O 

c. 

-75.58596 

-76.01075 

-76.19685 

A(42  —  43) 

• 

-40.5 

— 

•61.7 

— 

-68.3 
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TABLE  30.  TOTAL  (fa  ktrtmt)  AND  RELATIVE  (fa  kcml/aMl)  ENERGIES  OF  INTERMEDIATES 
AND  TRANSITION  STRUCTURES  FOR  THE  CYCLIZATION  OF  1-NITROPROPENE. 


Stractwaa 

Sym. 

HF/3-21G 

ibi. 

nt. 

HF/6-31G* 

aba. 

id. 

MP2/6-31G* 
aba.-  rel. 

CjHjNOj 

30 

c. 

-318.71225 

0.0 

-320.54169 

0.0 

-321.48544 

0.0 

30  (30  — »  37) 

c. 

-318.64118 

44.6 

-320.44925 

58.0 

-321.41821 

42.2 

37 

c. 

-318.71557 

-2.1 

-320.51239 

18.4 

-321.45489 

19.2 

44  (37  -  46) 

c. 

-318.68262 

18.6 

-320.46408 

48.7 

-321.40995 

47.4 

45 

c. 

-318.72281 

-6.6 

-320.51789 

14.9 

-321.46039 

15.7 

40  (46  — »  47) 

Cl 

-318.66661 

28.6 

-320.44728 

59.2 

-321.40630 

49.7 

47 

Cl 

-318.75295  - 

-25.5 

-320.53116 

6.6 

-321.45667 

18.0 

48 

c. 

-243.21070 

-244.58850 

-245.33411 

H]0 

c. 

-75.58596 

-76.01075 

-76.19685 

A(42  —  43) 


-27.4 


-42.7 


-46.6 
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SUMMARY  AND  CONCLUSIONS 

In  this  project  we  have  demonstrated  the  viability  of  nitronic  acids  in  die  smallest 
alkanes,  alkenes,  and  aromatics. 

The  patent  nitro  «*  ori-nitro,  keto  enol,  and  imine  enamine  tautomeric  processes 
were  investigated  at  various  levels  of  theory,  including  Gl.  The  energy  difference  between 
die  tautomeric  pairs  reduces  significantly  when  large  basis  sets  are  used  in  conjunction  with 
extensive  inclusion  of  electron  correlation  effects. 

Tantomerism  is  enhanced  in  die  conjugated  systems  which  were  calculated  at  MP4/6- 
31G*.  These  systems  allow  for  intramolecular  [1,5]-H  transfers.  Hydrogen  bonding  is 
present  in  the  2-nitrovinyl  alcohol  and  amine  and  their  nitronic  acids  have  trans  O4-N-O-H 
conformations.  In  contrast,  that  of  1-nitropropene  has  a  cis  form.  The  energy  differences  of 
die  nitroethylenes  with  die  ds  nitronic  acids  are  much  larger.  Cyclizadon  of  1-nitropropene 
can  occur  via  its  aci-nitro  tautomer  and  is  exothermic  after  dehydration.  Cyclizadon  of 
nitrovinyl  amine  (via  a  trans  nitronic  add)  is  less  likely  and  not  possible  for  the  alcohoL  In 
these  cases  decomposition  is  more  likely  as  the  N-OH  bond  is  relatively  weak. 

Tantomerism  in  orAo-substituted  nitrobenzenes  is  a  [1.5J-H  transfer  process  that  leads 
exclusively  to  nitronic  adds  with  a  ds  (N-N-O-H  conformation.  These  systems  woe 
calculated  at  MP2/6-31G*.  The  energy  differences  between  the  nitro  and  oci-nitro  tautomers 
are  significant  with  the  smallest  difference  for  o-nitrophenol  and  the  largest  for  o-nitrotoluene. 
The  nitronic  adds  have  sizable  kinetic  stabilities. 

Tantomerism  in  CHS-,  NHj-,  and  (XI-  substituted  2 ,6-dinitrobenzerves,  calculated  at 
HF/6-31G*.  is  similar  to  that  in  die  mononitro  derivatives  with  subtle  differences.  The 
buttressing  effect  of  the  two  nitro  groups  enhances  tantomerism  in  tire  phenol  and  toluene 
derivatives,  but  hampers  that  in  die  aniline  derivative  possibly  because  of  the  double  H- 
bonding  of  the  amine  group. 

Tantomerism  in  picric  add,  picramide,  and  TNT  was  calculated  at  HF/6-31G*  at 
which  level  the  nitronic  adds  have  significant  kinetic  stabilities.  With  the  additional 
polarizing  4-nitro  group  die  energy  difference  between  the  tautomers  (and  barriers)  is  reduced 
further  for  the  phenol  and  toluene  derivatives  but  increases  for  the  aniline  analogue. 
Tantomerism  is  by  far  the  easiest  for  picric  add,  (which  is  experimentally  also  the  most 
sensitive  to  impact  of  the  three),  with  an  energy  difference  of  25  kcal/mol  and  a  barrier  of  34 
kcal/mol.  The  nitronic  adds  of  picramide  and  TNT  are  significantly  less  stable,  Le.  48  and  39 
kcal/mol,  respectively,  but  follow  the  experimental  sensitivity  order  of  these  compounds. 

MP2  correction  on  the  [1,5]-H  transfer  barrier  for  tautomerism  of  TNT  gives  a  value  of  52 
kcal/mol  which  is  similar  to  the  experimental  induction  activation  for  thermal  decomposition. 

These  studies  demonstrate,  conform  to  experiment,  that  picric  add  by  means  of  its 
nitronic  add  is  more  sensitive  than  picramide  and  TNT.  The  theoretical  energies  on  the  latter 
compare  with  experimental  data.  The  studies  also  show  that  high  levels  of  theory  are  needed 
for  accurate  calculations.  They  further  establish  that  competing  decompositions  may  exist  for 
TNT  and  must  exist  for  picramide. 
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